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Abstract 

Based upon the lattice Dirac operator satisfying the Ginsparg- Wilson relation, we investigate 
canonical formulation of massless fermion on the spatial lattice. For free fermion system 
exact chiral symmetry can be implemented without species doubling. In the presence of 
gauge couplings the chiral symmetry is violated. We show that the divergence of the axial 
vector current is related to the chiral anomaly in the classical continuum limit. 



PACS: 11.15.-q, ll.15.Ha, 12.38.Gc 



1 Introduction 



The discovery of lattice Dirac operators satisfying the Ginsparg- Wilson (GW) relation |2l 
0] enables us to implement exact chiral symmetry on the lattice ^ E] without suffering from 
species doubling [H]. The overlap Dirac operator [3 , for instance, not only possesses all the 
desired classical properties such as the continuum limit and locality [7] but also reproduces 
correct chiral anomaly [Hlini- The lattice Dirac action with the exact chiral symmetry may 
be considered as a correct starting point of the nonperturbative studies of gauge theories 
with massless fermions lllj. 

In lattice gauge theories it is usual to employ euclidean path integral formalism. The 
main concern there is to compute various physical quantities nonpertubatively. To investigate 
formal aspects of the theory concerning states and operators it is suitable to work in canonical 
approach. To the author's knowledge, canonical treatment of GW fermion has only been 
investigated by Creutz, Horvath and Neuberger ^2] • They found an interesting characteristic 
structure of the spectra of energy and axial charge leading to axial anomaly. In this paper 
we will pursue the validity of their approach further and establish the chiral anomaly in the 
classical continuum limit. 

On the full eucliean lattice, it is possible to give action with exact chiral invariance 
at the classical level The fermion measure, however, is not chiral invariant but gives 
raise to nontrivial Jacobian fHl E] • It reproduces chiral anomaly in the classical continuum 
limit [HI El- In the canonical approach time is a continuous coordinate and only the spatial 
coordinates are discretized. This causes a problem in constructing chirally symmetric action. 

To put this more precisely, we consider a Dirac operator* 



where Vq is the time component of the ordinary covariant derivative in the continuum and 
Dd~i is the spatial part of D. It is assumed to satisfy the GW relation. If we choose Dd-i to 
be the overlap operator, species doublers can be avoided. The full Dirac operator D, however, 
does not satisfy the GW relation even in the absence of gauge couplings. Nevertheless it 
is possible to define exact chiral symmetry in the free field case. The chiral transformation 
there is analogous to the one given in Ref. The axial charge of a free GW fermion 

depends on momentum. In the physical momentum region it is almost constant as in the 
continuum theory. The separation of positive and negative axial charges decreases as the 
energy increases. It vanishes at the doubler momenta. In other words the states with positive 
axial charge and those with negative one are smoothly connected with one another at the 
boundaries of the first Brillouin zone [T^ . 

In the presence of gauge couplings, the chiral transformation becomes field dependent. It 
is no more a symmetry of the action. This is quite different with the conventional euclidean 
approach. The violation of the chiral symmetry is not a bad news. We expect it from the 
beginning since the chiral symmetry must be broken in any regulalized theory, otherwise we 
could not reproduce the chiral anomaly. We show that certain gauge fields induce asymmetric 
flow in the spectrum that leads to nonconservation of the axial charge. This is the origin of 
chiral anomaly. 

The interpretation of axial anomaly that external gauge fields can generate asymmetric 
level shifts of the negative energy states in the Dirac sea is well-known in the continuum 

*See Sect. Olfor notation. 



D = i^^Vo + Dd-i , 
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theory Jl]^ Similar analysis has been be carried out also for Wilson fermion. In this case 
the chiral symmetry, however, is broken break by the Wilson term even in the absence of 
gauge couplings. The axial anomaly for Wilson fermions was obtained in Ref. and its 
physical picture was investigated in Ref. |l5j . 

This paper is organized as follows. In the next section we consider free GW fermion and 
introduce exact chiral sjTiimetry. In Sect. El we examine the effects of couplings with the 
gauge fields on the chiral symmetry and give the axial current divergence. In Sect. |3]we 
describe the canonical approach explicitly for a two-dimensional system and examine the 
conservation of the axial charge for an adiabatically changing external electric field. In Sect. 
El we compute the classical continuum limit of the axial current divergence in an arbitrary 
smooth background and show that the chiral anomaly is correctly reproduced. Sect. El is 
devoted to summary and discussion. We collect notation in Appendix |X1 The technical 
detail in carrying out the momentum integrations to compute the anomaly coefficients is 
given in Appendix iBl 



2 Exact Chiral Symmetry 

In this section we consider canonical formulation of a lattice Dirac theory in an arbitrary 
even dimensional space-time with continuous time t and d — 1 discretized lattice coordinates 
X. The signature is assumed to be minkowskian and the spatial part is hypercubic regular 
lattice with lattice spacing a. The massless fermion action can be generically written as 



Sf = a' 



/ dtiJ(^^J'^^+D,^^^^P , (2.1) 



where D^-i stands for the spatial part of the lattice Dirac operator. Unlike euclidean for- 
mulation ip and ip are related by the Dirac conjugate 

^ = . (2.2) 

The conventions for the metric and 7-matrices are summarized in Appendix ^ 

To keep the fermion massless while avoiding species doubling, we employ the following 
construction analogous to the overlap operator 

/^d-i = -f l + 7d+i-^ 1 , (2.3) 



where Hd-i is the hermitian Wilson-Dirac operator given by 



Hd-, = Id+iiaD^w'^ - m) , D^^'^^ = J^i^'di - rdt) . (2.4) 



fc=i 



The (9f and are, respectively, symmetric and antisymmetric difference operators defined 
by ()A.4jl and is nothing but the Wilson-Dirac operator with the difference along the 

t There are several unpublished works. See references cited in jl5l Quantized field theory approach 
was investigated in Ref. [Tz] . 
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d-th (euclidean time) axis omitted. This choice is consistent with the cubic symmetry of the 
underlying lattice. The r and m are parameters. They must be so chosen that the theory 
contains no doublers. Here we assume r > 0. 

The Dirac operator ()2.3j) satisfies the following relations 



-fd+iDd-i + -Dd-i7d+i = aDd-ijd+iDd-i , 7°-Dd-i + Da-ij^ = aDa-ij^Dd-i , (2.5) 

where the first is the Ginsparg- Wilson relations as in the case of overlap Dirac operator 
on full euclidean lattices. The second one is specific to the present canonical approach. It 
can be seen by noting the fact that D^-i commutes with j'^'-fd+i = i^'^'j^ . . .7*^"^. It also 
satisfies the following hermiticity 

Dl_,-f' = ^^Dd-i , D\_,^d+i = Id+iDd^i . (2.6) 

From the action ()2.1|) we define single-particle hamiltonian h by 

h = —f'Dd^i . (2.7) 

The hermiticity of h can be seen from the 7" hermiticity ()2.6|) . 

In the absence of the interaction with gauge field the action ()2.1|) is invariant under the 
lattice chiral transformation 

dip = iejd+i ^1 - '^^d-ij ^ , = ietjj ^1 - ^-^d-i^ Id+i , (2.8) 

where e is an arbitrary real parameter. These chiral transformations are consistent with the 
Dirac conjugate (j2.2j) and the hermiticity of Dd~i- 

The invariance of the action under the chiral symmetry also implies that the axial charge 
operator qd+i defined by 

qd+i = Id+i - \Dd-^ . (2.9) 
commutes with h. Furthermore, h and qd+i satisfy ^2] 



"\^ + g^^, = l. (2.10) 



4 

These can be checked directly by using ()2.5|1 and ()2.(j|l . We thus see that (a£'/2, Qd+i) with 
E and Qd+i, respectively, the eigenvalues of h and qd+i can be regarded as a point on the 
unit circle and unlike the continuum theory the axial charge depends on the energy of the 
state. 

To see the doubling problem we go over to momentum representation 

d-l 

ij{t,x) = M(p)e*(P"-^*) . (px = ^pV) (2.11) 

k=l 
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To parametrize the eigenvalues of h and q^+i it is convenient to introduce d-dimensional 
orthonormal coordinates n by 



d-i 

— psinap^ , Till = — p^r ^^(cosap^ " 1) + i^ 



k=l 
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A sin^ ap*^ + 1 r ^^(cosap''' — 1) + m 1 . (2-12) 
\ k=l 



k=l 

These define a mapping from T'^~^, the 1st Brillouin zone —-< p^ < - with the opposite 
boundaries identified, to 5"'"^. The energy eigenvalues E and the axial charge Qd+i are given 
by 



|£ = ±yi±!^, Q,,,=±^i^. (2.13) 

Doublers appear at = or vr/a except for the origin p'^ = if = — 1 is simultaneously 
satisfied. For < m/r < 2 the energy E vanishes only at the origin . 

In quantum theory we impose the equal-time anticommutation relations 

{Mt,x),i^l{t,y)} = 6i,p6,,y/a''-' , x), y)} = x), V'Jlt, Z/)} = ,(2.14) 

where a and (3 are spinor indices. The hamiltonian, the fermion number and the axial charge 
are given by 

n = a^-^J2^^h^ ^ (2-15) 

X 

Q = a'^-i^^t^^ (2.16) 

X 

Qd+1 = a^-i^V^^g^+i^. (2.17) 

X 

By expanding the field operators in terms of plane wave solutions with definite energy and 
axial charge, we can introduce creation and annihilation operators. In Sect. lU we will 
explicitly carry out this in (1 + l)-dimensions. 



3 Coupling with Gauge Fields 

We now introduce the coupling with gauge field. This can be achieved by simply replacing 
the time-derivative with the covariant derivative Vq = d/dt + Aq and the differences dk, 
with the covariantized operators and defined by (jA.5jl . The fermion action ()2.1|1 is 
then replaced by 

Sp = a'^-i J2 fdt ^{ijoVo + Dd-M ■ (3.1) 
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This is invariant under the gauge transformation 

ipitjx) K{t,x)ip{t,x) 

Aait.x) k{t,x)AQ{t,x)k\t,x) - dtK{t,x)k\t,x) (3.2) 
Uk{t,x) A{t,x)Uk{t,x)A'^{t,x + ak) , 

where k is the unit vector in the fc-the direction and Uk{t,x) is the hnk variable associated 
with the hnk (x, x + ak). 

If we consider the gauge fields as dynamical, we must also introduce their kinetic part 
to the action. The magnetic part of the field strengths can be defined from the standard 
plaquette variables 

Pfc;(t, x) = Ukit, x)Ui{t, X + ak)Ul{t, X + ai)Uj{t, x) . (3.3) 

The electric field is given by 

Ek{t, x) = -VkAoit, x) + -Uk{t, x)Ul{t, x) . (3.4) 

CI 

The Wilson action for the gauge field takes the form 

X ^ X k<l ^ 

where g is the coupling constant. 

In general the overlap Dirac operator becomes singular for gauge fields where the hermi- 
tian Wilson-Dirac operator has a zero-mode. This implies that the fermion action ()3.1|) is 
not a well-defined functional over the entire space of the link variables. It is well-known in 
the euclidean path integral formulation that the overlap Dirac operator is well-defined for 
gauge field satisfying so-called admissibility condition 

\\l-Pki\\<e, (3.6) 

where 1 1 ■ ■ ■ 1 1 stands for operator norm and e is some positive constant depending on the 
parameters r, m and the dimensionality d jTj. By excising nonadmissible configurations 
the space of lattice gauge fields acquires nontrivial topological structure. The admissibility 
condition helps to define the chiral and gauge anomalies precisely in the full euclidean lattice 
theory. 

Practically, nonadmissible gauge fields can be avoided by modifying the gauge field action 
so that they are decoupled from the system [TDj. Incidentally, in (1 + 1) dimensions there 
is a parameter region of r and m where the overlap Dirac operator Di can be defined for 
any gauge field. This case is of great interest since the system corresponds to the lattice 
regularization of the massless Schwinger model. In the remainder of this paper we shall 
consider the lattice gauge fields as background external field. 

The chiral transformation ()2.8|) . however, is not consistent with the coupling of the gauge 
field. The spatial part of the fermion action ()3.1|) is invariant by construction. However, 
the term containing time-derivative violates the symmetry. The difficulty comes from the 
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noncommutativity of Vq and D^-i- In fact under the variation ()2.8|) with e being an arbitrary 
local parameter the action changes as 



+ [e,D,_i]7d+i(-l + |^d-i + ^(^7o2^o + ^d~i)) h ■ (3.7) 



If e is a constant parameter, the first and the last term in the integrand are absent. Hence 
the chiral symmetry is violated if [Dq, Dd-i] is nonvanishing. Incidentally, we can retain the 
chiral invariance for static gauge fields in the temporal gauge Aq = 0. 

The variation ()3.7|) of the action leads to the axial current divergence relation 



where we have introduced the following notation to simplify the expressions 

Dd-i^ = Dd-Mt, x) , ^Dd-i = ^Dd^iit, x) . (3.9) 

The axial charge density appearing in the time derivative of this expression coincides with 
()2.17|) . The second term on the Ihs of ()3.8|) corresponds to spatial divergence of the axial 
current. In fact the integral of this term over the spatial lattice vanishes. This can be seen 
by noting the relation like 



What we have seen is that in the presence of the couplings with an external gauge 
fields the commutativity of Dd-i with the time derivative is lost and consequently the chiral 
symmetry ()2.8p is violated. In Sect. El we will show that the violation of the chiral symmetry, 
the rhs of ()3.8|) . which is naively of 0(a), reproduces the chiral anomaly in the classical 
continuum limit. 

Since the chiral transformation ()2.8j) is not an exact symmetry, we may consider naive 
chiral transformations 

Sip = iejd+itp , H = ieip'jd+i, (3.10) 

which is not a symmetry of the action as well. We close this section with a comment on 
what happens if one uses naive chiral transformations instead of ()2.8|1 in the evaluation of 
chiral anomaly. The variation of the action 1)3.111 under the naive chiral transformation gives 
the following axial current divergence 

^{^7%+i^} - #(^d-i - Dd-i)i' = -ia^jj D d^i-fd+iDd-i^. (3.11) 

Compared with ()3.8|) . we arrive at an apparently different form of chiral anomaly if we employ 
fl3.10|) as the chiral transformation. These two, however, differ only by a total divergence of 
some gauge invarinat current k^^ as can be seen from 

- la^Dd-ild+iDd-i^ = -^ii^old+i[Vo, Dd-i]^ + + d^k' (3.12) 



(3.^ 
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where k'^ is defined by the relations 

= |V^7o7d+i^d-i^, d*ik' = --^'^(Dd_i-Da-i)ld+iDd-i'^ (/ = 1,2,3). (3.13) 

The axial currents corresponding to the transformations ()2.8|1 and ()3.10|) differ only by a 
gauge invariant current of 0(a) and give essentially the same chiral anomaly in the classical 
continuum limit. 



4 Ginsparg- Wilson Fermion in (1 + l)-dimensions 

In this section we apply the formalism developed in the preceding sections to the case of 
(1 + l)-dimensions. We consider a finite periodic lattice of size L = a{2N + 1), where we 
have assumed that the number of sites is odd. This is only a technical assumption to make 
the arguments simple. 

In momentum space the eigenspinors satisfy the Dirac equation 

-f'Eu{p) + Di{p)uip) = , Di{p) = -(7^ni + + 1) , (4.1) 

a 

where Di{p) is the overlap operator in momentum space. In (1 + l)-dimensions the eigen- 
spinors can be more conveniently parametrized by a pseudo-momentum p defined by 

sin ap 

sm ap = —Til = 



^sin^ ap + (r(cos ap — 1) + m)^ 



r(cosap — 1) + m 
cosap = —Ti2 = . (4.2) 

y sin ap + {r{cosap — 1) + my 

For < m/r < 2 these define one-to-one mapping from —it /a < p < n/ato —it /a < p < tt / a. 
We consider p as continuous function of p beyond the first Brillouin zone. In particular we 
have 

~27r 27r , , 

p+ — =p+ — . 4.3 
a a 

One easily see from ()2.12|) that the choice r = m = 1 leads top = p. In (l + l)-dimensions this 
parameter choice is special in the sense that the overlap operator Di becomes identically 
equal to the Wilson operator. This can be verified directly by establishing the relation 
Hi = 1. Remarkably, it holds true even in the presence of gauge couplings. 

The eigenspinor for E ^ is automatically eigenspinor for the chiral operator defined 

2 ap ap 
by ()2.9j) . We thus find plane wave solutions for E = Ep = - sin -— and Qs = Xp = cos ■ 



a 2 ^ 2 



URp{t, x) = — w^(p)e^(P-^^*) , ur{p) 



( ap \ 

cos — 
4 ^ 
ap 

V— T/ 



(4.4) 
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Figure 1: The energy (•) and the axial charge (o) are plotted as functions of momentum for 
r = 0.8 and m = 1.2 with N = 5 (left). They are conveniently described as points on a unit 
circle (right) [12]. Diametrically opposite points correspond to a pair of eigenspinor uji{p) 
and ul{p). 



and for E 



2 ap 
-Ep = — sm — and Q3 
a 2 



cos ■ 



ap 



ULp{t,x) 



UL{p)e 



i{px+Ept) 



ul{p) 



( . ap\ 
sm — 

ap 



(4.5) 



The dispersion relations E = ±Ep are easily recognized as the lattice analog of E = ±p 
known for the continuum chiral theory. Unlike the fermion number the conserved axial 
charge depends on the momentum. For the physical region \p\ « n/a the axial charge is 
Q3 = ±Xp ^ ±1. It approaches to zero at the boundaries of the 1st Brillouin zone \p\ ~ ir/a. 

The momentum p is usually restricted to lie in the 1st Brillouin zone. In the present 
chiral theory the period of the spectra of energy and axial charge is not 27r/a but is ATv/a. 
Furthermore, the wave functions ()4.4j) and ()4.5p up to an overall sign are transformed into 
each other by the translation p —>■ p + 2'n'/a. As will be clear shortly, it turns out to be 
convenient to double the Brillouin zone and use a single wave function defined by 



Up{t, x) 



u{p)e 



i{px—Ept) 



u{p) 



( ap 
COS — 

4 

. ap 



(4.6) 



By noting the relation ()4.3|1 we see that up to an overall sign Up is periodic in p with a period 



Air /a and Up = Urp for — vr < ap < t: and Up 



U 



Lp—2-K /a 



for vr < < Svr. The relation 



among the momentum, the energy and the axial charge are shown in Fig. ^ 
The free Dirac field can be expanded in terms of Up as 

V'(t,a;) = ^ apUp{t,x) , 

0<ap<A-K 



(4.7) 
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where ap is a creation or an annihilation operators satisfying the anti-commutation relations 

{ttp, a^,} = 6p.pf , {ap, ttp'} = {al, a^,} = . {0 < p, p' < An /a) (4.8) 

We see from ()4.6|) that aj, creates a particle of energy Ep ioi Q < p < 2tx / a and ap create 
an anti-particles of energy \Ep\ = —Bp for 2TT/a < p < An /a, where we have assigned the 
zero-energy states as particles for simplicity. 

The hamiltonian H, the fermion number Q and the axial charge Q3 are given by 

H = a^i!^hip= ^ Epa^ap , (4.9) 

X 0<ap<47r 
X 0<ap<47r 

Q3 = a^i'^'q^ip = 5Z ^p4^p • (4- 

X 0<ap<47r 

These are not normal ordered with respect to the creation and annihilation operators and 
the energy and the fermion number of the Dirac vacuum are nonvanishing. However, the 
axial charge of the Dirac vacuum vanishes because of the cancellation of the axial charge 
between the occupied negative energy states. 

We now introduce U{1) lattice gauge field {AQ(t,x),Ui(t,x)) in (1 + l)-dimensions and 
investigate the behavior of the axial charge. In the remainder of this section we assume that 
the gauge potential is real. The electric field defined by ()3.4|) is given by 

£(t,x) = -diAo(t,x) - -dtUUt, x)Ui(t, x)^ , (4.12) 

a 



which is invariant under the gauge transformation ()3.2j) . In this section we are only interested 
in a spatially uniform electric field generated by the gauge field 

Aoit,x) = , Ui{t,x) = e-'"^W . (413) 
Then the electric field is given by 

£{t) = -A{t) . 

Since the gauge field is translation invariant in the spatial direction, we can solve the Dirac 
equation in momentum space. In the presence of the background gauge field the free hamil- 
tonian — 7°Di(p) is modified to — 7°Di(p — A{t)). Then the Dirac equation in momentum 
space is simply given by 

= -i'D,{p-A{t)Mt,p) , (4.14) 
At this point, we further assume that the A{t) changes adiabatically from A{ti) = to 

where is an integer satisfying |ng| < A^. We take t/ — ti sufficiently large so that the 
adiabatic change of A(i) is possible. In this situation the transitions between ub, and are 
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suppressed and a state that is initially in an eigenstate of the hamiltonian keeps staying in 
the corresponding eigenstate of the hamiltonian at an arbitrary time. Thus the solution to 
the Dirac equation ()4.14|) satisfying (f{ti,p) = u{p) can be solved in terms of the eigenspinor 
of -7°L'i(p- A(t)) as 



'^(t,p) = u{p — A(t)) exp 
where u{p) is given by ()4.6|) and Ep 



2 ap 
— sm — 
a 2 

electric field the free Dirac field ()4.7p evolves to 



Vl 



apu{p — A(t)) exp 



0<ap<47r 



Ep_A(T)dT , (4.16) 
Under the influence of the background 

(4.17) 



px- Ep_AiT)dT + 6p 



The phase dp is chosen so that ip{t, x) reduces to the free field (j4.7j) for t <ti 
For t >tf the Dirac field approaches to 



ijj{t,x) 



iqx 



Vl 



^ apu{p - q)e 



i((p~q)x-Ep-gt+5'p) 



(4.18) 



0<ap<47r 



where Sp is some irrelevant phase. 

The overall factor e^'^^ is periodic on the lattice if the condition ()4.15|) is satisfied. It can 
be removed by carrying out a time-independent gauge transformation 



ij{t,x) 



-iqx 



ip{t,x) 



(4.19) 



This transforms iplt, x) {t > tf) into a free field. As an external gauge field, it is not necessary 
to impose the condition ()4.15j) and any q = A{tf) is allowed. In general e*^"^', however, is 
not periodic on the lattice and cannot be regarded as a gauge transformaiton on the lattice. 
Only those satisfying ()4.15|) can be eliminated by the gauge transformaiton ()4.19|) . 

What we found is that the state with initial momentum p is carried over to the state 
with momentum p — q due to the electric field exerted on the system. It is well-known in 
continuum theory that this kind of spectral flow results in the violation of the conservation 
of the axial charge. To see this we compute the vacuum expectation value of the axial charge 
aX t = tf. It is given by 



(0|Q3(t/)|0) 



Y: cos^^(0|aX|0) 

0<ap<47r 




COS ■ 



2 

a{p + q) 



cos ■ 



a{p - q) 



^ cos 

TT—2aq<ap<TT 



a{p + q) 



cos ■ 



a{p - q) 



(4.20) 



Q<ap<2aq 



where |0) is the Dirac vacuum satisfying ap|0) = for < ap < 27r and use has been made 
of ()4.3|1 . In the last equality we have assumed g > 0. The contributions from the modes 
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lying near the boundary of the Brillouin zone (vr — 2aq < ap < tt) are canceled out and only 
the modes around the origin {0 < ap < 2aq) contribute to the violation of the axial charge 
conservation. We thus find 

(o|Q3(t/)|o) = - ^°^y- (4-21) 

—aq<ap<aq 

The g < case can be analyzed in a similar way. 

To relate ()4.2H) with the axial anomaly we consider the classical continuum limit a ^ 
with L = a{2N + 1) kept fixed. For the modes contributing to the nonvanishing axial charge 
ap 

we have cos — 1 as a — > 0. We thus find 

(0|Q3(t/)|0) = -2ng = --a(2Ar+l)g=-aV [ ' S{t)dt , 



where use has been made of £(t) = —A(t). Since (0|Q3(tj)|0) = 0, we finally obtain 

/■*/ . 1 r^f 

/ dt{0\Q3{t)\0) = 2— dx dtS{t) . (4.22) 
Ju ""^^ 27rJo Ju 

This coincides with the integrated axial anomaly relation. 

The nonconservation of the axial charge implies the violation of the chiral symmetry. In 
the next section we investigate the chiral symmetry in the presence of an arbitrary back- 
ground gauge field and see how the chiral anomaly is reproduced in the classical continuum 
limit. 



5 Chiral Anomaly 

In the presence of a background gauge field the chiral symmetry is violated. We will compute 
the violation of the axial charge by examining the classical continuum limit of the vacuum 
expectation value of the operator appearing in the rhs of ()3.8|) . It is a local function of the 
gauge field and is denoted by 2A{t, x) 

A{t,x) = -^{ilj{t,x)'yo'yd+i[Do,Dd-i]4j{t,x)) 

= -Uim lim tr7rf+i7°[Po,^d-i](T^(t,a;)V^(t',a;')) , 

4 t'^t x'^x 

where tr stands for trace over spin and internal indices and T denotes time-ordering. The 
two-point function {Tip(t,x)ip{t',x')) satisfies 

(^7°I?o + Da-i){T^{t,x)^{t',x')) = i6{t - t')'5x,.'/«'"' • (5.1) 
Hence A can be written as 
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where the ie prescription is made exphcit. For simphcity we assume that the size of the 
lattice is infinite. The rhs of this expression can be simphfied by making use of the relation 



(5.3) 



After a bit of algebra, we find 

A{t,x) = 



4Qd-2 

i 



dE 


r d^-^p 




l-n (27r)<^-i 


dE 


r d'^-^p 








1 



+ ie 



tr 7d+i7 



T^o, Id- 



+1 



Hd-i 



HI, 



X- 



(5.4) 



where we have introduced the notation O — ^-w^jaQ^rpxia f^j. operator O. In the last 
equality we have used a rescaling trick aE — > E. 

To carry out a systematic expansion by the lattice constant a we further assume that 
74o(t,a;) approaches to a smooth function as a — >■ and the link variable are given by a 
smooth gauge potential A}^{x) as 



[/fc(t, x) = V exp 



Then the expansions of Vf and are given by 



a / dsAkit^x \ a{\ — s)k) 



(5.5) 



aVf = i sinp'' + a cosp^Vk + O(a^) , 
aVj^ = cosp'' - 1 + sinp'^Pfc + O(a^) , 



(5.6) 
(5.7) 



where T>k = dk+A^ is the covariant derivative in the continuum theory. For later convenience 
we use n^, nif^ (/x = 1, • • • , d) defined by 



-psmp^ , -p{r^{cosp'' - 1) +rnj , 

k 



rrik = pcosp 



md = p 



(5.8) 



(p ^ = , sin^p^ + (r y~^(cosp^ ~ 1) + ^) ? A; = 1, 
V fc k 



■ ,d-l 



In terms of and the hermitian Wilson-Dirac operator if^-i and its square can be 
expanded as 



p-fd+iHd-i ^-f-n + Hd + Vi , p^Hj_, = 1 + V2 + V3, 



(5.9) 



13 



— rn) — 7 ■ n — 



where 7 ■ = 'j^rik and V^s are given by 

A;=l 

k k 

= ia ^(7''mfc + rnk)Vk + 0{o?) , 

k 

V, = -aV5^(Vff + /(ar5^V^ + m)'-l 

k k 

= 2ia^{mk + rndinkVk + 0{a^) , (5.10) 

k 

V, = .raV5^7'=[Vf,Vf] + iaV5^7V[Vf,Vf] 

k,l k,l 

k,l 

In the last hne use has been made of the field strength Fki = [Dkj'Di]. The reader might 
think that we should expand V2 to order a^. In the actual computation the order term in 
V2 does not contribute to A at least in two and four dimensions. 
From ()5.9|) we obtain 



7d+i#d-i/y^^ii = h-n + nd + Vi)/^/l + V2 + Vs . (5.11) 
If we define V, which is of order a, by 



V = za7°Po + id+iHd-i/ \l Hj_^ - ^ ■ n - , (5.12) 
we can expand ()5.4p as 



d-2 



A{t,x) = ^^('=) + 0(a) , (5.13) 

fc=0 

where A^'^^ is given by 

■ r+co 7 rp fir Ad—^rr. 

= l,m(-l)'+'j^ If j_^A-^,,-,,,,^''lV,„V](SVfS . (5.14) 

The S is the free propagator 

S = -— ^ . (5.15) 

7OE + 7 ■ n + 1 + rirf + ze ^ ^ 

It is now straightforward to compute A^''\ 
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5.1 (1 + l)-dimensions 

The computation of ()5.4|) is rather simple in two dimensions. From the general argument 
given above only the lowest order term A^^^ survives in the limit a ^ 0. We thus obtain 



. . N i f^^ dE r dp n 
A(t,x) = / — / — tr737° 



Vo,V^-\{l-n + n2)V2 



7°^ + 7- n + l + n2 + «e ^ ' 

-^cse'^'^tr F^, + 0(a) , (5.16) 
47r 



where C2 is a numerical coefficient given by 



This integral can be evaluated by a change of variable z = psinp as explained in Appendix 
IBI We find from ()B.9|) C2 = 1 for < m < 2r and C2 = for m < 0, m > 2r. This implies 
that the chiral symmetry breaking ()5.1fjj) approaches to the chiral anomaly of the continuum 
theory for < m < 2r. Note also that the result is consistent with ()4.22|1 . 

5.2 (1 + 3)-dimensions 

In four dimensions we must compute A^'^^ for k = 0,1,2. Inserting ()5.12j) into the rhs of 
()5.14|1 and carrying out the trace over the 7-matrices and the energy integral, we obtain 

2 



p ^ J j,k,l 



n (27r)3 ^2(1 + n4) 



^(2) ^ 1 r 1 




16 (27r 



We do not reproduce the computation here. Instead, we mention some properties useful in 
getting some feeling about the results. Firstly, all the contributions which are potentially 
diverging in the limit a —>■ vanish by virtue of the 75 in the trace. Similarly, the terms of 
order a^ in Vi, V2 and ia'-f^Vo in V, which should be considered in the systematic expansion 
of V in the lattice constant, can be ignored. Then we have only to consider the terms 
explicitly given in ()5.10p . Secondly, the trace over the spinor indices followed by the energy 
integral yields momentum integrals of functions of m^, and 714. Typically they take the 
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forms 

j,k,l \ p J j,k,l 

r d^pJ2^'''rn,m,minfO'^,,{n,) = f [ J] m J ^ e^'^'O;,,^) , (5.19) 

-'-^ j,k,l -f--^ \ p / jAl 

-'^'^ j,k,l,p -'^'^ \ p / \ p / j,k,l 

In moving from the Ihs to rhs of these expressions we have used the fact that n4 is symmetric 
in the momentum variables. 

What remains to show is the evaluation of the momentum integrals. This is done in 
Appendix El Using ()B.13j) . we find that the chiral symmetry breaking is given by 

At,x) = -^e^-^P'^tiF^^Fp^ + Oia) . (5.20) 

Since n+ = 1 for < m < 2r, this completely agrees with the chiral anomaly in the 
continuum theory. We thus establish that the canonical description of the Ginsparg- Wilson 
fermion correctly reproduces the chiral anomaly in the classical continuum limit. 

6 Summary and Discussion 

We have investigated canonical formulation of massless Dirac theory on the spatial lattice. 
In the absence of gauge couplings the theory possesses an exact chiral symmetry of the 
Liischer type while avoiding species doubling. The axial charge operator commutes with 
the hamiltonian and the conserved axial charge of a particle depends on its momentum. In 
the classical continuum limit the momentum dependence of the axial charge disappears and 
the states with opposite chirality are decoupled with each other. On the lattice, however, 
the spectra of energy and axial charge are smooth periodic functions of momentum with 
a period Air/a and the gap between positive and negative axial charges disappears at the 
corner points of the first Brillouin zone. Then transitions bewteen states with opposite signs 
of axial charge may occur by the gauge couplings. This is responsible for the axial anomaly 
as was noted in Ref. [12]. 

In the presence of gauge couplings the chiral transformation depends on the gauge fields. 
The violation of the chiral symmetry in our canonical approach is simply due to the fact that 
the axial charge operator does not commute with the hamiltonian. Our computations show 
that by taking account of the fermion loop effect the Ward-Takahashi identity for the broken 
axial charge conservation correctly reproduce the well-known anomalous conservation law in 
the classical continuum limit. 

One might think that the GW fermion with the chiral transformation ()2.8|) looses su- 
periority to the Wilson fermion with the naive chiral transformation ()3.10p in the presence 
of gauge couplings. The naive chiral symmetry, however, is violated even in the absence of 
gauge couplings whether one uses the Wilson fermion or the GW fermion. In the case of the 
modified chiral trasnformation it is only broken at the gauge couplings. We expect that the 
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breaking of the modified chiral symmetry is more controllable than that of the naive one. 
This is the virtue of using the GW fermion. 

The interpretation of the axial anomaly as the spectral flow for the adiabatic change of the 
gauge field could be extended to higher dimensions as in the continuum theory ^1 HHl ■ 
In four dimension we can consider a uniform external magnetic field. If the gauge field is 
time independent, we can still define conserved axial charge. The energy spectrum can be 
parametrized by the momentum along the magnetic field. We expect that the dispersion 
relation similar to the two-dimensional case discussed in Sect. |3] arises per flux quantum for 
sufficiently smooth gauge field Applying uniform electric field parallel to the magnetic 
field would induce axial charge nonconservation proportional to the total magnetic flux. 

In the euclidean path integral formulation the chiral anomaly is directly related to the 
index of the lattice Dirac operator and gives a topological invariant when summed over 
the lattice. Topologically nontrivial structure of the configuration space of lattice gauge 
fields emerges by imposing admissibility condition. This happens even in two dimensions. 
In the canonical approach the axial charge changes continuously. The configuration space 
of the lattice gauge field is divided into sectors by applying the condition that the system 
approaches, up to gauge transformations, to a free field asymptotically. 

It is natural to ask whether the formalism can be applied to chiral gauge theories like the 
standard model. To achieve this it is necessary to define chiral fermions. In the euclidean 
path integral approach fermion field and the conjugate field are independent variables and 
chiral fermions can be defined by using different chiral projection operators for fermion field 
and the conjugate field. In the canonical approach this does not work even for the free theory 
since the conjugate field is related to the fermion field by the Dirac conjugate. This difficulty 
could be circumvented by carefully choosing the fermion contents so that the would-be gauge 
anomalies be canceled. 

Finally the issue of quantizing the gauge field lies beyond the scope of the present paper. 
It is interesting to pursue the understanding of nonperturbative aspects of QCD such as the 
^^-vacuum in the present approach. 

We thank Yoshio Kikukawa and Hiroshi Suzuki for valuable discussions. This work is 
supported in part by the Grant-in- Aid for Scientific Research from the Ministry of Education, 
Cuhure, Sports, Science and Technology (No. 13640258, No. 13135203). 

A Notation 

The metric is assumed to be rj'^^ = diag(l, — 1, ■ ■ ■ , — 1). The 7-matrices satisfy 

{7^7'^} = 2r/^^ 7°^ = 7°, 7^^ = -7', 

7.+i = ^^-Vy---7'-\ 7^1 = 1. (A.l) 
In (1 + 1) -dimensions we employ 

, _ _i _ / 1 1 ,1 _ „ „2 _ / ~1 \ _ ,,0,1 _ 3 



7" = = [I oj , 7^ = = [I J , 73 = 7%^ = a' . (A.2) 

Using the forward and backward difference operators dk and given by 

1 . 1 

dkip{t, x) = -{ip{t, X + ak) — il){t, x)) , dlipit, x) = -{ipit, x) — ipit, x — ak)) , (A. 3) 
a a 
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we define the symmetric and antisymmetric difference operators by 



dl = l{dk + dl), dt = ^id,-dl), (A.4) 

where k stands for the unit vector along the k-th spatial coordinate axis. Note that 9f ~ 
0(aO) and ~ 0(a). 

Covariant difference operators are defined by 

Vk'ipit, x) = -{Uk(t, x)i'(t, X + ak) - tp(t, x)) , 
a 

Vlipit, x) = -{ip{t, x) - Uk{t, X - akyip{t, X - ak)) . (A. 5) 

a 

The symmetric and antisjTiimetric covariant differences are defined by 

Vf = ^(Vfc + VD, ^t = \{^k-Vl). (A.6) 

B Momentum Integrals 

The momentum integrals similar to ()5.17p and ()5.18|) already appeared in the evaluations of 
chiral anomalies in full euclidean lattice theories. They can be carried out by considering a 
continuous map from p'^ G T'^^^ to G S'^^^ as in Ref. [S], where the coordinates on S'^^^ 
are defined by 

= -Uk = psinp^ , z'^ = -rid . (/c = 1, ■ ■ ■ , d - 1) (B.l) 

The volume form d'''~^z of the {d — l)-disk D'^'^ defined by {z^Y + " " " + {z'^^^Y < 1 is related 
to the volume form d'^'^p of the momentum space by 

d'^-h = d'^-^pJ{p) , (B.2) 

where J{p) is the Jacobian 
Here we are interested in the integral 



p \ p ^ 



+ I± = f d''-'p9iTnd)find)[]Jmp^(nd-rJ2^] , (B.4) 



p \ p ^ 



where 6 is the unit step function. Let be the domain in the momentum space 

that is mapped into the upper (lower) hemisphere of S'^^^ with z'^ > {z'^ < 0). Then the 
insertion of ^(=F^d) effectively restricts the momentum region to D±. Since = =FVT^^r^ 
with r = a/ (z^)"^ + ■ ■ ■ + {z^~^)'^ on D±, we obtain 



=F/i± / d z 



d~l 



18 



where fi± is defined by 

/i± = sgnJ(p) . (B.6) 

{p\zd{p)=±l} 

If there is no point on (D^) that is mapped to the north pole with z'^ = 1 (the south 
pole with z'^ = —1), then /x+ = = 0). 

The computation of fi± is similar to the counting of the doubler modes in Ref. We 
consider the points (p^, • • • G T'^~^ with being either or tt. They are mapped to 

the north pole {z'^ = 1) or the south pole {z'^ = —1) of S'^~^ by (jB.lll . For the parameters r 
and m satisfying 2n < m/r < 2(n + 1) (n = 0, ■ ■ ■ , d — 1), the points {p\, ■ ■ ■ ,pi~^) where 
at most n entries are tt are mapped to the north pole, otherwise to the south pole. We thus 
find 



s=0 ^ ^ s=n+l ^ ^ 



We see that fi± satisfy 

fi^ + ^_ = . (B.8) 

If we assume ^ ) ~ ^ ^'^^ s<OoTS>d— 1, (jB.7j) becomes valid for any integer n. It 

yields /i± = for m < or m > 2(i. 

In two dimensions we find fi± = ±1 for < m/r < 2 and /i± = for m < or m/r > 2. 

Carrying out the remaining radial integral for / = y — -^^ — - and noting ()B.4|1 . we obtain 



■'I 



, l + n2 ( \ —4 for < m/r < 2 _ , 

V 2 V mi/ 10 for m < 0, m/r > 2 



To find the anomaly coefficients appearing in ()5.18|1 we need only two types of momentum 
integrals /i and I2 defined by 



^2(1 + 714) W J\ p^P 



rrip 
p ^ 



In four dimensions ^± are given by 

{0 for m < 0, m/r > 6 
1 for < m/r < 2, 4 < m/r < 6 (B.12) 
-2 for 2 < m/r < 4 

Using ()B.4jl and ()B.5j) and carrying out the radial integrations, we find 

IGtt IGtt 

^1 = — :^^^+ ' ^2 = -r^^^+ ■ [b.u) 

6 15 
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